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The sensitivity of the Rabi oscillations of a resonantly driven spin-| system to a weak and slow 
modulation of the static longitudinal magnetic field, -Bo, is studied theoretically. We establish the 
mapping of a weakly driven two-level system with modulation onto a strongly driven system without 
modulation. The mapping suggests that different regimes of spin dynamics, known for a strongly 
driven system, can be realized under common experimental conditions of weak driving (driving 
field Bi <^ Bg) upon proper choice of the domains of modulation frequency, aj,„, and amplitude, 
B2. Fast modulation cj„, ^ SIh, where is the Rabi frequency, emulates the regime of driving 
frequency much bigger than the resonant frequency. Strong modulation, B2 3> Bi, emulates the 
regime Bi ^ Bq. Resonant modulation, Um ~ ^r, gives rise to an envelope of the Rabi oscillations. 
The shape of this envelope is highly sensitive to the detuning of the driving frequency from the 
resonance. Theoretical predictions for different domains of B2 and tj„ were tested experimentally 
using NMR of protons in water, where, without modulation, the pattern of Rabi oscillations could be 
observed over many periods. We present experimental results which reproduce the three predicted 
modulation regimes, and agree with theory quantitatively. 

PACS numbers: 42.50.Md,76.20.+q,76.60.-k 



I. INTRODUCTION 

Ever since the oscillations of the population of Zeeman 
levels in a resonant ac magnetic field were predicted the- 
oretically by Rabff', they have been experimentally ob- 
served in numerous media and in various spectral ranges 
from the radio to the optical. The revival of interest in 
Rabi oscillations during the past decade has been fueled 
by the fact that they can now be measured in a single 
two-level system (qubit). In fact, the observatiorpHlSl of 
high-quality Rabi oscillations in certain isolated two-level 
systems such as an exciton in quantum dotPSl^ a Joseph- 
son junctiorPt^, a single spin in a dolP^^^!^ or a vacancy 
gpjjJi3H15l jg viewed as evidence that the corresponding 
qubit can be coherently manipulated. 

Conventional experimental conditions for observation 
of the Rabi oscillations are: 

(i) the driving frequency, w, is in resonance with the split- 
ting, Az = g^oBg, of the spin levels. Here g is the g- 
factor and ^0 is the Bohr magneton. 

(ii) in-plane driving field, xi?iCoswi, is much smaller 
than the static longitudinal field, zBg. This condition 
guarantees that the Rabi frequency, ^Ir = g^oB^, is much 
smaller than A^. Under the conditions (i), (ii) the clas- 
sical Rabi resullfi] 
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FIG. f: Schematic illustration of magnetic resonance in a 
modulated longitudinal field. The Zeeman splitting, Az, os- 
cillates in time with a small magnitude, £„, and frequency. 



which was empty at t = 0, applies. Here 



(2) 



is the oscillation frequency, and (5 = Az — wisa detuning 
of the drivin g freq uency from resonance. 

R is knowrpSEIthat almost sinusoidal Rabi oscillations 
can also be excited in the domains of parameters where 
the conditions (i) and (ii) are not met. One exampld^ 
is the multiphoton Rabi oscillations which develop when 
the driving frequency is close to 



,(p) - 
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for the oscillating population of the upper Zeeman level, where p is positive integer. Another exampl^^ is realized 
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in the domain of frequencies 



(4) 



Existence of the domain Eq. Q requires a strong driv- 
ing field Bi ^ Bo, which is quite unusual for magnetic 
resonance. 

The main point of the present paper is that the en- 
tire spectrum of regimes of the Rabi oscillations can be 
realized when conditions (i) and (ii) still apply, but the 
longitudinal field contains a small slowly oscillating com- 
ponent, zB2{t), where B2 ^ i?o- In other words, we 
will assume that the Zeeman splitting, A^, has a time 
dependent correction. 



s{t) = gfioB^it), 



(5) 



and study how the classical result Eq. ([!]) is affected 
by this correction. We will demonstrate that even when 
the magnitude of the correction is much smaller than 
and it oscillates with a frequency much smaller than A^, 
see Fig. [l] its effect on the Rabi oscillations can still be 
dramatic. To prove this statement in Sect. II we map a 
weakly and resonantly driven two-level system with mod- 
ulation onto a strongly driven system without modulation. 
In sections III, IV, and V we make use of this mapping to 
study how different regimes of spin dynamics emerge in 
different limits of modulation, namely, fast (compared to 
flu) modulation, strong (compared to fl^) modulation, 
and near-resonant modulation (with frequency « ^r)- In 
Sect. VI we report the results of the experimental test of 
theoretical predictions. NMR measurements on protons 
in water have advantage that, without modulation, the 
pattern of Rabi oscillations could be observed over many 
periods. By applying the modulation pulse with vary- 
ing frequency and magnitude, we were able to reproduce 
the above three regimes of spin dynamics. Concluding 
remarks are presented in Sect. VII. 



II. MAPPING ONTO A STRONGLY DRIVEN 
SYSTEM WITHOUT MODULATION 



equation by introducing, instead of C_^_l and the 
combinations 



C, 1 +C_i 



c_ 



(8) 



The system of equations for the new amplitudes reads 



jA+ = —^A^ + coscji. 



iA^ 



'-Aj^ — rtjiA^ cosujt. 



(9) 



(10) 



Upon expressing A- from Eq. ^ and substituting it 
into Eq. (10), we arrive at the sought equation 



A, 



—ibjfla sin ujt + cos^ uit 



A+=0. (11) 



Assume now that the modulation is present, but, due 
to the weakness of the ac field, fin ^ A^, the system 
Eqs. ^, can be treated within the rotating-wave 
approximation (RWA) . This amounts to the replacement 
of coswt by ie-*"* in Eq. and by ^e*"* in Eq. 
Then, upon standard transformation into the rotating 
system 



these equations take the form 



iD^ 



iD_ 
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(12) 



(13) 



(14) 



As a next step, we express D_i from Eq. ( 13 ) and substi- 
tute it into Eq. (14). This leads to the following second- 



order differential equation for D. 



The time evolution of the amplitudes C_|_i and C_ 1 of 
two spin orientations is governed by the system of equa- 
tions 



iC.i 



iC_ 



A2 +£(<) 



C, 1 -I- f2HC_i coswt. 



(6) 



-e{t) 



C_i + ri/;C_|_i cosLut. (7) 



In the absence of modulation, e{t) = 0, it is convenient 
to reduce this system to a single second-order differential 



Mt) , (s + sr + nl 



D^i = 0. 



(15) 



At this point we assume that the modulation is sinu- 
soidal. 



e{t) ~ e„, cos w,„t, 



(16) 



and make the key observation that, for S = 0, Eq. (15) 
reduces to Eq. (11) upon replacement 



A, 



2n, 



UJ. 



(17) 
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This mapping provides an important insight into the ef- 
fect of modulation of the longitudinal field on magnetic 
resonance. Indeed, Eq. (11) captures the time evolution 
of the populations |C_|_i(tJf^ and |C_i(i)p when the ac 
drive is not weak, i.e., > A^, so that RWA does 
not apply. On the other hand, Eq. (15 1 describes the 



spin dynamics with weak ac drive, but in the presence 
of the modulation. Solutions of Eq. (Ill, the struc- 



ture of which is dictated by the Floquet theorem, were 
analyzed in a grea t number of papers starting from pio- 
neering work a^^ l I Therefore, the approaches developed 
for a driven spin-i system without modulation can be 
utilized for the case when the modulation is present. 

For example, in the absence of modulation, resonant 
drive of the spin-i system corresponds to the condition 
w = Az. Then from Eq. (17 1 we conclude that the 



effect of modulation on the Rabi oscillations is most pro- 
nounced when the modulation frequency, w„„ is close to 
the Rabi frequency, O^,. This sensitivity was previously 
pointed out in Refs. 



A less obvious consequence of the mapping Eq. (17) is 
that, upon increasing the modulation strength, e„, the 
Rabi oscillations become sensitive to the modulation for 
modulation frequencies close to 



,(p) - 



2p+l' 



(18) 



Condition Eq. ( 18 ) implies that the modulation affects 



the Rabi oscillations when the modulation period con- 
tains an odd number of the Rabi periods. This condition 
is an analog of the condition for multiphoton resonances 
in a driven two-level systeml^which take place at ac driv- 
ing frequencies uj ~ uj^p\ where oj^^^ is given by Eq. (jsj. 
Despite the fact that Eqs. (11) and (15) can be for- 



mally mapped onto each other, the physical phenomena 
that they describe are vastly different. For the most 
prominent example, w ~ Az, Eq. (11) describes the 



Rabi oscillations of the populations of the Zeeman levels, 
while Eq. ( 15 ) describes slow envelope of these oscil- 



lations, which emerges due to modulation, e(t), of Az- 
Formally, the difference stems from the initial conditions. 



which must be imposed on the solutions of Eqs. ( 11 ) and 



( 15 ). If, for example, at i = the spin points down, then 
the initial conditions for Eq. (11) read 



A+(0) = 1, i+(0)=z 



Az 
2 



(19) 



For the same initial state, the initial conditions for Eq. 
( 15 ) have the form 



Z?+i(0) = 0, Z),i(0) 



(20) 



Below we study the spin dynamics, P_|_i(<), in different 
domains of modulation frequencies and magnitudes. 



III. FAST MODULATION: aj„ > Qr. 

Qualitatively, it is clear that weak modulation e„ ^ 
w„ averages out if the modulation frequency is much big- 
ger than ilfl. Then the Rabi oscillations remain unaf- 
fected. Nontrivial modification of the Rabi oscillations 
takes place when both w„ and e„ are much bigger than 
ilfi. For conventional magnetic resonance this regime cor- 
responds to w » Az and fl^ ^ Az, see Eq. (17), i.e.. 



Bi ^ B„, which is exotic. Fast modulation, on the other 
hand, is fully compatible with Bi <Si B„. 

For simplicity we will consider the case of zero detun- 
ing, 6 = 0. To analyze the limit 3> ^Ir it is convenient 
to introduce the new variables 



2 2 ^ ' 

^+i-i^+i(i)e^ 



- sin LOjrit 



(21) 



and to rewrite the system Eqs. (13), (14) as 



^R 1 



^R 



- sin cu™ t 



D 



(22) 



If we now replace exp (±-^j™- sinw^t) by the time average, 
' '^here Jo{z) is a zero-order Bessel function, the 
system Eq. (22) will readily yield 



D_i =cos ( ' ^ -is\-a {^J'i^ ' * 

(23) 

(24) 



so that 



2 / ^R T / £,1 



P+i(0=sin^(^Jc 



The latter expression suggests that the effect of fast mod- 
ulation, w„ 3> ^R, on the Rabi oscillations is the reduc- 
tion of their frequency by a factor Jo{^f^ ■ The reduction 

is significant when the modulation amplitude, e,„, is of 
the order of cj„. 

The remaining task is to demonstrate that, for ar- 
bitrary modulation strength, e„, the condition w„ 3> 
VI R justifies the replacement of exp (±^j™- sin w„t) by 

"^o(ff") ■ ^^^^ purpose we consider the correction to 
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FIG. 2: Fast modulation regime, a;„ ^ Q^. Modulation 
causes slowing down of the Rabi oscillations, which are plot- 
ted with dashed lines from Eq. ( |24[ ) for tj„, = 5Q,r and three 
values of the ratio — : (a) 0.5, (b) 1.5, (c) 2. Correction 



Eq. ( 25 1 gives rise to the fast component of the Rabi oscil- 



lations. The probability, P^i{t), with fast component taken 
into account is shown with full lines. 



D_i_ 1 coming from higher harmonics of exp (^j^ sinw^i) . 
It is convenient to cast this correction in the form 



2W„ ,yo 



2 



D^i + i sin ( -^Jol —— I t 



- sin c^^jje. 



where we introduced a small parameter 



cos 770, (25) 



(26) 



There is a small prefactor in front of the integral. Still 
one has to check that the integral does not grow at large 
LOrnt- Suppose that w„t = 2-kN . Then the integration 
over N modulation periods can be reduced to a single 
integral from to 27r in which cos 7](j) is replaced by 



cos /sin27r?7iV 
cos ri[(f> — 2'Kn) = — ^ — ( — — — - — f-cos 2'KriN + 1 

n=0 



tan TTTj 



sin 7y0 / 1 — cos 2'ki]N 
tan TTrj 



sin 2TTriN 



(27) 



We see that for N < i]^^ the first term grows with TV. 
Moreover, this term is much bigger than the second term. 
However, the correction Eq. (251 is determined by the 
second term in Eq. (27). The reason is that the first 
term does not depend on (p with accuracy rj. On the 
other hand, the part of the first term which is indepen- 
dent of (f) does not contribute to the correction Eq. (25 1. 
This is because the integral from this part is zero. The 



second term in Eq. (27 1, by virtue of smallness of rj, can 



be replaced by ^ sin nrjN, so it remains finite at large 
iV. This proves that, even for strong modulation, the 
correction Eq. (25) is ^n/uj^. In terms of the proba- 
bility P^i(i), the correction Eq. (25) gives rise to a fast 
component, as illustrated in Fig. [2j 



IV. STRONG MODULATION: e„ > > 

For strong modulation, 3> the term in Eq. 
(15) is small compared to which suggests that Rabi 



oscillations do not develop. However, if the modulation is 
slow enough, and the detuning, (S, is much smaller than 
rjfl, the term will be dominant during short (com- 
pared to ijJ^^) time intervals near 



(28) 



when e{t) passes through zero. Within these intervals we 
can set t = + ti and expand e{t) with respect to ti. 
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Then Eq. ( 15 ) assumes the form 



D+i = 0, 



(29) 



where we assumed for definitiveness that k is odd. 

Solution of Eq. (29) can be expressed in terms of 
parabohc cylinder functions^'^ 



where the characteristic time, r, is given by 

1 



and parameter v is defined as 

v = -i —. 

4e„ uj„. 



(30) 



(32) 




FIG. 3: Strong modulation regime, 2> fij? 3> o^m- Spin dy- 
namics persists only during short intervals \ t—tk\ r <^ i^It'^'' 
(31) around the moments, t^, when modulation passes through 



zero. Probability, P^i, is plotted from Eq. (36 1 versus di 
mensionless time {t ~ tk)/T for three values of parameter 
Eq. (32 1, describing the strength of the ac drive. Green, red, 
and blue curves correspond to the values \u\ = 0.1, \v\ = 0.3, 
and \v\ — 0.7, respectively. 



Expansion of e(i) is justified under the condition, w,„t ^ 
1, which reduces to e„ ^ With regard to parameter, 

I/, it can be presented as a product \iy\ — (^^'^^ (^)' 
The first factor of this product is small, while the second 
factor is large. This means that in the domain e,„ 3> 
^ u!^ the value \v\ can be both large and smaU. The 
magnitude of determines the angle of the spin rotation 
as e{t) passes through zero. To see this, assume that at 
moment t = t^. the spin points down. Then the constants 



/i) /2 in Eq. ( 30 ) can be found from the initial conditions 



Eq. (20), yielding 



215^0) 



(33) 

where 'DKO) is the derivative of 'D^(z) at z = 0. The 
dependence of f ^0) on v can be established from the in- 
tegral representation of the parabohc cyhnder functio'^ 



■DUO) = 



1/2 



2Z./2 



sm 



TTiy 



(34) 



Here T{z) is the gamma- function. The probabihty 
f+i(ti) contains |I?[,(0)p, which, using the properties 
of the gamma- function, can be simplified to 



PUO)|' = |^^|sinh 



7r|i/| 



(35) 



Substituting Eq. (35)intoEq. (33 ), we arrive at the final 
resuh 



1 



4sinh(^ 



(36) 



In Fig. [31 the dependence P+i(ti) is plotted from 
Eq. (36) iOT different values of the parameter \v\. We 
see that, starting from zero, P^i{ti) oscillates with 
ti, and, finally, saturates at ^ t at some finite 
value, P_|_i(oo). To find P^i(cx)) we use the large-z 

m 



asymptote of the parabolic cylinder functiorP^, 'D^(z) '-^ 
z'^exp and the fact that the asymptotes of T>^{z) 

and I?^(— z) differ by exp(«7ri/) — exp(— 7r|z/|). This al- 
lows to find the saturation level of the numerator to be 
4sinh^^^^^ exp ^— Subsequently, the saturation 
level of Pj_i assumes the form 



1 



-7r|i/| 



(37) 



The meaning of the result Eq. (37) is transparent. If 
the driving field is weak, so that ^ 1, we have 
P_i_ 1 (oo) 1, which implies that the spin pointing down 
at t = tk retains its orientation throughout the entire in- 
terval t^^i — = If the driving field is strong, so 
that I J/ 1 ^ 1, the spin, during the short time <C — , will 



rotate to the position P_^_i ss P_ 



will point 
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FIG. 4: Schematic illustration of the second-order Rabi oscil- 
lations emerging as a result of slow modulation of the static 
field, -Bo- Driving resonant ac field, splits the quasienergies 
by Qr; modulation with frequency cj„ « leads to an ad- 
ditional splitting. 



along the a;-axis, and spend in this position the remaining 



part of the interval 



V. WEAK RESONANT MODULATION: 
A. Qualitative picture 



In Sect. II, based on the mapping Eq. (17), we con- 



cluded that Rabi oscillations are strongly sensitive to the 
modulation of the longitudinal field when cj„ is close to 
0„. It is illustrative to discuss the effect of modulation 
near this condition using the language of quasienergies. 
Within this language, without modulation, the presence 
of a resonant field with frequency w « results in a de- 
generacy of quasinenergies corresponding to Zeeman-split 
levels, see Fig. |4j Conventional Rabi oscillations emerge 
as a result of lifting this degeneracy (avoided crossing); 
the magnitude of the splitting of quasienergies is Qr. 

With periodic modulation, the Hamiltonian in the 
quasienergy representation remains time-dependent and 
the corresponding eigenstates are characterized by 
second-order quasienergies. When the frequency of mod- 
ulation is close to ^Ir, the first-order Rabi-split quasiener- 
gies become degenerate, leading to their additional split- 
ting, see Fig. |4j The magnitude of this additional split- 
ting is controlled by the modulation amplitude. In the 
time domain, this additional splitting manifests itself via 
second-order Rabi oscillations in the form of an envelope 
of the primary oscillations. 

As we show in the next subsection, the shape of this 
envelope is very sensitive to the detuning, a;„ — Hr, of 



modulation frequency from the Rabi frequency, and to 
the detuning, S, of the driving frequency from the Zeeman 
splitting, Az. 



1.0 
0.8 

0.6 
'0.4 
0.2 
0.0 



20 40 60 80 100 

(a) Ofi t 



(b) 



1.0 
0.8 
0.6 
0.4 
0.2 
0.0 

1.0 


0.6 
' 0.4 
0.2 
0.0 




20 40 60 80 100 




20 



40 



60 



(C) 



80 100 



FIG. 5: Evolution of the Rabi oscillation pattern with de- 
tuning of the modulation frequency from the Rabi frequency. 
Spin dynamics, P^i{t), under the resonant drive, 5 = 0, is 



plotted from Eq. (44 1 for different dimensionless detunings. 



k: 0.5 (a), 1.0 (b), and 1.5 (c). The amplitude of modulation 
is £,„ =0.1 Qr. Upon increasing k the envelope develops (a), 
reaches maximum (b), and gradually vanishes (c). 
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B. The shape of the envelope 



where we introduced a dimensionless deviation 



Near the condition a;,„ = 17^ Eq. ( 15 ) can be solved in 



the resonant approximation, which amounts to keeping 
only two terms in the Floquet expansiorP^: 



i\t 



do + d_ 



(38) 



where A is the Floquet exponent. Substituting Eq. (38) 
into Eq. ( 15 ) and equating the resonant terms, we get the 
following system of coupled equations for the amplitudes 

do, d^i 



ni 



S' + ^-4y] do = -{5- uj^)e^d_ 



splitting of the Floquet exponents 



2(u;,„ - a;o)a;o 



(41) 



of the modulation frequency from oscillation frequency, 
LUo and took into account that \uj^ — Wo| <C Wq. The 
steps leading from Eq. pS] ) to Eq. (401 are the same 
as the steps leading from Eq. (Ill to the conventional 
Rabi oscillations Eq. ([l]). The difference arises when one 
substitutes A+, A_ into the system Eq. (39), finds two 



ni + S' + ^ -4{\-u;„y d_, = -{5 + uj^)e^do. (39) 



For vanishing modulation, e„ 0, both left-hand sides 
in the system Eq. (39) turn to zero at the degener- 
acy point A = ^a;„, when the condition: a;„ = = 
^/6^~+~^}^ is met. At finite the system yields the 



linearly independent solutions, D^i{t) and D i{t), and 

. . . ^ . 

requires that their sum satisfies the initial conditions Eq. 
(20). After that, upon calculating the occupation 



2 n- o ^ o 



(42) 



± 



4 Wo 



one arrives at the following generalization of Eq. ([T]) to 
(40) the case of modulated longitudinal field 



1 



(5(Kfi„ -I- S) 



4 Wo 



+ 6) . 2 /t^o I / , , -, \ \ 

— K i H , sm K — V + 1 \t 

! ) 2uji{^i^^^~k) \2 4wo V ^ 

» (^ _j^ sin^ i^-^^'-^U^ y^)) t 
2w2(y^^^ + K) V2 4wo V Jj 



It 



(43) 



We start the analysis of Eq. ( 43 ) from the case of zero 



detuning, 6 = 0, when it simplifies considerably 



P,i{t) = 

+ 2 ^ ^ k2 + 1 



sm I — + -^K]t 



2 1 ^R , £r, , 

sm + 



2(Vk2 + 1 + k) V 2 4 V J J 



(44) 



This expression is consistent with Ref. [22]where a related 
quantity, namely the absorption of an external field by a 



resonantly driven spin-i system in a modulated longitu- 
dinal field, has been studied. 



The last two terms in Eq. (44) describe the second 



order Rabi oscillations since their frequencies are split 
in agreement with qualitative picture Fig. |4j Naturally, 
taking the limit — >■ 0, we recover the conventional 

Rabi oscillations, P_^_i (t) = sin^ (^^^) j for any k. A non- 
trivial feature of Eq. (44 1 is that exactly at a;„, = iln, 
when K = 0, the last two terms cancel each other, i.e., the 
Rabi oscillations are unaffected by the modulation. Grad- 
ual development of the second-order oscillations upon in- 
creasing K is illustrated in Fig. [5] We see that the effect 
of modulation of longitudinal field is most pronounced at 
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dimensionless detuning k « 1, where the modulation of 
the Rabi oscillations is complete, see Fig. [5]d. Upon fur- 
ther increasing k, the effect of modulation vanishes above 
K = 4. 

Beyond the resonant approximation'^, the difference. 



u!^ — u!o, in Eq. (41) acquires a correction ^ 
which is an analog of the Bloch-Siegert shiflP^. 

Solving Eq. ( 15 ) in the resonant approximation is also 
permitted when a;„ is close to cj^-' = -^fj- Then the 
two resonating terms in the Floquet expansion are dg 
and d_^2p+l)■ The splitting of the Floquet exponents, at 
the degeneracy point, A — 2(2p+i) ' takes place in the 
(2p + l)-th order in parameter e,„/riR. The magnitude 
of the splitting can be obtained by readjusting notations 
in the expression for the width of multiphoton resonance 
obtained in Ref. [T7| 



Ap) 



In the domain 




(45) 



UJ„ 



) ' Eq. ( 44 ) remains valid 



upon the replacement e,„ and describes second- 

order multiphoton Rabi oscillations. Naturally, if the 
modulation is not purely sinusoidal but contains a har- 
monics, say I, the envelope of the Rabi oscillations near 
a;,„ — fla/l will already develop in the lowest order. 

In the absence of modulation, finite detuning, S, of the 
driving frequency from the resonant frequency, Az, leads 
to the reduction of the period of the Rabi oscillations 
and to the reduction of their amplitude, see Eqs. ([T]), 
([2]). Remarkably, in the presence of modulation, near 
the condition a;„, = ^Ir, finite d causes a new qualitative 
feature in the shape of the Rabi oscillations. Moreover, 
the effect of finite 6 depends strongly on its sign. 

In Fig. [6] we plot from Eq. ( 43 1 the evolution of the 
Rabi-oscillation pattern as S increases in the positive di- 
rection, while parameter k is kept constant. We see that, 
while for small enough S, the prime effect of modulation 
is still a slow envelope of the oscillations, upon further 
increasing of S, the envelope gradually transforms into 
the oscillations of the baseline. At large enough 6, Fig. 
[6];, these oscillations of the baseline dominate the Rabi- 
oscillation pattern. 

As 6 increases in the negative direction, see Fig. [?) the 
pattern of the Rabi oscillations evolves differently. At 
small |i5| modulation gives rise to the envelope, similar 
to the case of small positive S. However, at a certain 
\S\, see Fig ^p, the magnitude of the oscillations does 
not change with time at all, suggesting that the effects of 



two detunings, uj„ 



and, 



Az, cancel each other. 



Only upon further increase oi \5\ the magnitude of oscil- 
lations decreases significantly, and they proceed around 
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FIG. 6: Regime of weak resonant modulation. Evolution of 
the Rabi oscillation pattern with detuning, 5, of the driving 
frequency, lu, from the level spli tting, Az. Spin dynamics, 
P+i{t), is plotted from Eq. (43 1 for three positive values of 

0.3 (a), 0.6 (b), and 1.4 (c). The amplitude of modulation 
is Em = 0.15 Oh; dimensionless detuning of the modulation 
frequency from ila is k = 1.0. The pattern gradually evolves 
from envelope of the Rabi oscillations (a) to weak oscillations 
around the modulated baseline (c). 
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FIG. 7: Same as Fig. 6 for three negative values of —0.2 
(a), —2 (b), and —3 (c). The pattern evolves from envelope 
(a) to homogeneous oscillations (b). Weak modulation of the 
baseline emerges only at large detuning (c). 



the slowly oscillating baseline, as shown in Fig. [6]:. The 
plots in Figs. [6] [7]are drawn for the ratio e^/Qpt = 0.15 
and Em/^R — 0.2, respectively, however general features 
of the modulation pattern depend weakly on this ratio. 



VI. TESTING THEORY WITH NMR 
EXPERIMENT 

The introduction of the longitudinal modulation field 
eit) to investigate the general behavior of Rabi oscilla- 
tions is highly compatible with conventional experimen- 
tal NMR methods of inductive detection in thermally 
generated spin ensemble^^. For the large Bg field gen- 
erally required to achieve sufficient signal-to-noise ratio 
(SNR), the condition <^ Bg cannot be avoided, but 
the additional small-amplitude modulation field is eas- 
ily realized (Fig. [8^) . The three regimes elucidated in 
Sect. |ffl[ |IV[ and |V| were thus explored using straight- 
forward NMR of protons in water (gyromagnetic ratio 
1 = gfJ-N = 4.25775 kHz/G, where /iA? is the nuclear mag- 
neton). We needed only to take some additional care to 
ensure a stable and highly homogeneous rf driving field 
Bi across the sample, so that the pattern of Rabi oscil- 
lations could be observed over many periods. 

All experiments were performed in a horizontal-bore 2- 
Tesla superconducting magnet (Oxford Instruments). A 
conventional solenoidal single-coil transmit /recieve probe 
(5 turns, 1 cm diam and 2.5 cm long) was series-tuned 
with a capacitor to the proton Zeeman resonance at 
88.8 MHz. A 50-f2 resistor in series with these ele- 
ments provided a matching impedance to the transmit 
and receive amplifiers. This "low-Q" probe sacrifices 
some SNR for a robust and flat frequency response that 
is n egligi bly affected by the accompanying modulation 
fiel d^ I The water sample is centered in and occu- 




, I 

j Bq modulation 
j Helmholtz pair j 



FIG. 8: (Color online) a) A schematic of the NMR probe 
used in the experiments. A traditional NMR coil (B 1 j IS ac- 
companied by a Bo modulation Helmholtz pair that is coaxial 
with the -Bo field, b) Graphical description of the relation be- 
tween theoretical predictions and experimental NMR data. 
The theoretical predictions have the magnetization M pro- 
jected onto the z-axis of the Bloch sphere. The experimental 
data is the projection of the magnetization onto the xy-plane 
of the Bloch sphere. Eq. ( 46 1 relates the theory to experi- 
ment. 
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pies about 25% of the coil volume; it is contained in 
small tube made of PTFE (Tcflon'^'^), which minimizes 
the magnetic-susceptibility difference between the sam- 
ple and its immediate surroundings. Both the low filling 
factor and PTFE tube serve to enhance the _Bi -field ho- 
mogeneity across the sample. The modulation field was 
provided by a 5-cm-radius Helmholtz pair (coaxial with 
the main Bg field) wound on a form that has the probe 
coil at its center (see Fig. [8^). A NMR spectrometer 
(Tecmag Redstone, model HF2-1RX) with two indepen- 
dent transmission channels was used to transmit pulses 
to both the 88.8 MHz proton probe and the 0-100 kHz 
modulation coils, and subsequently to acquire and dig- 
itize the free-induction-decay (FID) signal generated in 
the probe coil. The rf pulse was amplified with a 
2 kW amplifier (Tomco model BTO2000-AlphaSA) con- 
ventionally designed for solid-state NMR, but whose high 
output power allowed the coherent nutation of the pro- 
ton spins through many Rabi-oscillation periods in a rel- 
atively short pulse time. The _Bo-modulation pulse was 
provided by a DC-50 kHz amplifier normally used for 
gradient coils in imaging applications (Techron AN7780) ; 
the inductance of the Helmholtz pair was matched to the 
gradient amplifier's output specifications in order to min- 
imize ring-down and associated cross-talk to the NMR 
coil. 

The FID signals were acquired on resonance (6 = 0) 
in single-shot fashion (no signal averaging), and the FID 
signal strength was plotted vs. B^ pulse length. Each 
data point displayed is the result of examining the corre- 
sponding FID to determine the magnitude of the trans- 
verse magnetization at a fixed time delay « 0.2 ms from 
the end of the B^ rf pulse. The B^ and modulation pulses 
are essentially applied simultaneously, with the modula- 
tion pulse nominally starting at e{t) — just as the B^ 
pulse starts. We assume that the i?i-pulse amplitude is 
fixed and the pulse length is linearly related to the nu- 
tation angle of the proton spins. Because of slight tran- 
sient changes in power delivered by the rf amplifier at 
the very beginning of the B^ pulse (minimized by using a 
low-Q probe), this linear relationship is most accurately 
observed for longer pulse times. Indeed, when we de- 
termine experimentally using non-modulated data, 
we rely on measuring the average period of the later- 
time oscillations. With this caveat in mind, the plots 
we show below represent Rabi oscillations in the three 
modulation regimes of interest. To decrease the overall 
data acquisition time, the water sample contained dis- 
solved copper sulfate (CUSO4) to reduce the longitudinal 
relaxation time to 100 ms. One must generally 
wait at least several times Ti between FID acquisitions 
to allow the magnetization to recover to its characteris- 
tic thermal-equilibrium value. Intrinsic transverse relax- 



ation (characterized by T2) is typically on the order of 
T-i in weakly interacting liquids. We observed transverse 
decoherence caused by residual inhomogeneity (charac- 
terized by T* < T2) in the B^ field, which manifests 
most clearly in our data as the overall decay of the non- 
modulated Rabi oscillations after many characteristic pe- 
riods (see Figs. 9pT|. 

We note that the predictions made in Sections |III||V| 
are formulated in terms of the projection of the magneti- 
zation onto the z-axis of the Bloch sphere (see Fig. [sja) . 
Rabi oscillations occur between the high- and low-energy 
states defined by B„ , where the low-energy state has mag- 
netization parallel to Bg. However, the observable in a 
conventional NMR experiment is the projection of the 
magnetization onto the a;?;-plane of the Bloch sphere. 
Noting also that the initial conditions at time t = for 
our predictions have P^i{0) = 0, whereas the experi- 
ments have P^i(O) = 1, a simple transformation of the 
z-axis prediction to the xy-plane can be accomplished by 



Pj_ ^ sin [arccos(2P^i - 1)] = 2. /P 



(46) 



In Figs. [9 |ll| we display the experimentally measured 
transverse magnetization data represented by P± and 
compare them to theory by transforming the predictions 
according to Eq. (46). The initial peak of the non- 



modulated Rabi-oscillation data, which appears in black 
in each figure, is used to define P± — 1 and to normalize 
the corresponding modulated data. 

We begin with the fast-modulation regime, uj^ 3> ^r, 
studied in Sect. |III[ The picture in the "rotating frame" 
that results from the application of the RWA, is that 
of an average nutation of the magnetization in the yz- 
plane, superimposed with much faster wiggles, which for 
small modulation amplitudes are transverse to the yz- 
plane. For larger modulation amplitudes, the wiggles 
move the magnetization appreciably along the surface of 
the Bloch sphere, giving rise to a time-average decrease 
in the component of the magnetization that is in the yz- 
plane and subject to a torque generated by B^. This 
leads directly to the slowing down of the Rabi nutation 
expressed by Eq. ( 24 p3. The wiggles also lead to the fast 
modulation component at uj^ , expressed in Eq. ( 25 ) , that 



is superimposed on the slowed-down Rabi oscillations. 
Fig. [9^ demonstrates that the corresponding experiments 
are sensitive to even the small decrease in the effective 
Rabi frequency ^IrJo (5™") that appears for weak mod- 
ulation. We were able to fit this data with Eq. ( 24 ) mul- 



tiplied by a a decaying exponential that accounts for the 
T* decay. The Rabi frequency and modulation am- 
plitude e„ were experimentally determined, respectively, 
by examining the non-modulated data and by measuring 
the current in the modulation coils. 
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FIG. 9: (Color online) Rabi oscillations in the fast- 
modulation regime (green) compared to non-modulated data 
(black). Transverse spin projection P± is plotted vs. pulse 
length; the latter is linearly related to the nutation angle, a) 
The modulation amplitude e„ is small compared with 
and the modulated data show the expected small upward 
shift in the oscillation frequency. A fit of the modulated 
data to Eq. ( 24 1 is shown by the dashed red line. The ex- 



perimentally determined values of the modulation frequency 
aj,„/27r = 41.30 ± 0.01 kHz, modulation amplitude £„/27r7 = 
4.96 ± 0.06 G, and the Rabi frequency nft/2TT = 8.26 ± 0.05 
kHz (determined from the non-modulated data) are input as 
fixed parameters. The fit includes multiplication of Eq. ( 24 1 



by 



to account for the overall decay; the fit then yields 



T2 = 1.28 ±0.02 ms. b) The modulation amplitude is large 
compared with ujm., and the oscillation frequency decreases 
significantly, in accordance with Eq. ( |24[ ). Here, uj ^/2tv = 
31.5 ± 0.01 kHz and r2H/27r = 3.10 ± 0.02 kHz are experimen- 
tally determined. A fit that incorporates the correction of 
Eq. (25 1 is shown in red, from which e„/27r7 = 13.0 ± 0.1 G 
is extracted; larger values of were difficult to measure ex- 
perimentally. The correction gives rise to the oscillations at 
ujm that ride on top of the much slower Rabi oscillations. The 
weaker Bi produces negligible T2 decay. 



Fig.|9}3 shows data in the regime of both fast and strong 
modulation £„„ ^ fl^, where corrections found in 
Eq. (25) become apparent. As predicted, the slowing- 
down effect on the Rabi oscillations now becomes pro- 
nounced, and the small modulation at frequency a;„ de- 



scribed by Eq. ( 25 1 rides on top of the envelope given by 
Eq. ( 24 1 . A peculiar effect seen in the data and predicted 
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FIG. 10: (Color online) Rabi oscillations in the strong- 
modulation regime (green) compared to non-modulated data 
(black). Transverse spin projection P± is plotted vs. pulse 
length; the latter is linearly related to the nutation angle. 
The lower part of the figure shows the complete data set. 
Data without modulation (black) yields a Rabi frequency 
r2R/27r = 17.7 ± 0.1 kHz. Data with modulation are shown in 
green; the applied modulation frequency was aj™/27r = 3.0±0.1 
kHz. The critical regions near e{t) = occur every 0.166 ms 
and are characterized by non-trivial oscillation behavior; the 
first such region is shown expanded in the upper part of the 
figure. These early-time data show the decreasing oscillation 
period characteristic of the parabolic cylinder func tion s. The 
results of a fit (red) of the early-time data to Eq. ( 36 1 yields 
e„/27r7 =7.0±0.05 G, which then leads to =0.9±0.1. An 
asymptote (dashed maroon) is calculated from Eq. ( |36[ ) and 
is also shown. Well away from the zero-crossings of e(t), the 
magnetization continues to nutate in a more regular fashion; 
that it nutates at all is the result of not sufficiently satisfying 
the limit ^ ^la. 



the nutation angle 7r/2, where the magnetization remains 
pinned on the xy-plane. In the fit of the data to Eq. ( 25 ) 
only e„ was used as a free parameter; for various tech- 



in the theory is the leveling-ofF of the magnetization near 



nical reasons, larger modulation amplitudes were more 
difficult to measure experimentally. The Rabi frequency 
rin. as determined from the unmodulated data was a fixed 
parameter in the fit. We note that the ^25% reduction 
(due to smaller B^) in the value of fij, determined from 
the non-modulated data appears to have eliminated the 
effects of T* decay in the unmodulated data of Fig. [9b 
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FIG. 11: (Color online) Rabi oscillations in the regime of weak resonant modulation, showing definite beats that depend on 
the parameter k. The Rabi frequency is found from the non-modulated data to be f2_R/27r = 16.7±0.7 kHz. The sensitivity of 
the fit to small changes in the modulation frequency and ex per imental limitations led us to float most of the experimental 
parameters and to perform two separate fits of these data to Eq. ( |44[ ) for the early-time (red dashed) and late-time (blue-dashed) 
portions of the overall acquisition. For the early-time fit, the extracted values are: 0^/2^ = 16. 5±0. 2 kHz, e„/27r7 = 0.91±0.04 
G, and aj,„/27r = 17.56±0.05, yielding k, = 0.54 ±0.14. for the late-time fit, the extracted values are: r2R/27r = 17.0±0.2 kHz, 
£„/27r7= 1.2 ± 0.1 G, and a;m/27r = 18.03±0.09, yielding k = 0.40 ± 0.12. The data exhibit the beating envelope characteristic 
of this modulation regime. 



as compared to similar data in Fig. [9^. 

In the strong-modulation regime e^^il^^ uj^ stud- 
ied in Sect. |IV[ the picture in the rotating frame is that 
the slowly sweeping modulation field brings the spins into 
resonance for only a short fraction of the modulation pe- 
riod; in the remaining time the modulation field is so 
large that the effective field in the rotating frame lies 
along the z-axis — the spins are essentially out of reso- 
nance, and do not nutate. Our experiments had e„ at 
most only a few times fij, and thus did not achieve the 
limit e,„ ^ to sufficient degree to turn off the nutation 
completely, even when \e{t)\ was near a maximum. How- 
ever, the seemingly complicated data seen in the lower 
half of Fig. [To] can still be understood to a significant de- 
gree. Non-trivial behavior of the Rabi oscillations occurs 
when e{t) is near zero with periodicity -k /(jJm ~ 0.166 ms, 
as per Eq. ( 28 1 . These critical regions of time are marked 



and the 
0.9 was calculated from Eq. (32). Using 



by dramatic changes in the oscillation behavior lasting 
about one Rabi period 27r/riji, whereas there is a rel- 
atively steady-state oscillation for the remainder of the 
modulation period. In general, the prediction given by 



Eq. (36 1 in Sect. IV can only be applied when the initial 



spin state is known, which is only true in our case for the 
critical region just after t = Q. 

Hence, the upper half of Fig. [TO] shows an expanded 
view of the first 0.05 ms of this data set, which was fit to 
Eq. ( 36 ) . The fit shown uses the experimentally deter- 
mined fifl = 17.2 ±0.1 kHz and determines e„ = 7.0 G as 



a free parameter. Using the values of f2;,, lo 
parameter \u 

this value the saturation level shown with dashed hori- 
zontal line was calculated from Eq. (36). We note with 
respect to the fit that e„ and the other parameters de- 
rived from it are extremely sensitive to whether the mag- 
netization reaches the a;?/-plane (Pi = 1). The precision 
quoted for e{t) may thus be somewhat underestimated. 
Nonetheless, the data and the fit clearly show the de- 
creasing oscillation period characteristic of the parabolic- 
cylinder functions on which the theory in this regime is 
based. 

The qualitative features of the data in Fig. [lO] away 
from the critical regions can also be understood. As the 
critical time period near a zero-crossing of e„ ends, the 
magnetization vector finds itself at some particular an- 
gle to the z-axis that would then essentially not change 
in the limit e„ 3> iln until the next zero-crossing. Even 
when this limit is not well satisfied, if the magnetization 
is close to being in the xy-plane, nutation about an effec- 
tive field that is mostly along the z-direction will produce 
only small second-order changes in the FID amplitude as 
the spins nutate: a clear example occurs from 0.2 ms to 
0.3 ms; here, the FID amplitude is very close to maxi- 
mum and there is a barely perceptible nutation. Larger 
oscillations occur for lower values of the FID amplitude 
when the magnetization is well away from the xy-plane, 
with the largest occurring between about 0.85 ms and 
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0.95 ms. We note that the dephasmg that occurs due to 
Bi inhomogeneity appears to be decreased by the strong- 
slow modulation. The modulation field has the eff'ect of 
continuously refocusing the spins in the rotating-frame 
yz-plane, analogous to the way a Hahn echo refocuses 
dephased transverse magnetization precessing about an 
inhomogeneous Bg. 

In our study of the weak-resonant modulation regime 
ila ~ u!m S> we focus on the limit where the B^- 
pulse is on-resonance with the Zeeman splitting, and the 
theoretical prediction for the spin dynamics is given by 
Eq. (44). In this regime we observe beats in the Rabi os- 



cillations; the depth of the modulation is determined by 



the parameter k. In Fig. 11 we show a typical example; 
a fit of the modulated data to Eq. (44) is qualitatively 



reasonable and clearly exhibits the beat envelope. Quan- 
titatively, a high-quality fit over the entire time interval 
is made difficult by the extreme sensitivity of k, which 
contains contains a small difference, (w„ — O^) <C to 
experimental errors in Because of the above compli- 
cations, we allowed most of the experimental parameters 
to become floating fit parameters and performed separate 
fits for two time intervals: "early" and "late", as shown 
in Fig. [TT] The values of k extracted from both fits were 
reasonably close to each other: 0.54 ± 0.14 for early and 
0.40 ± 0.12 for late fits. We attribute the discrepancy 
to the slow drift with time of the Rabi frequency and 
the modulation-pulse parameters. The evidence that this 
drift affects the fit can already be inferred from Fig. [9] for 
the fast-modulation regime. We note that, as was true 
in the strong-slow modulation regime, dephasing due to 
Bi inhomogeneity is suppressed by the modulation field. 



VII. CONCLUDING REMARKS 



is strong corresponds to e„ <C ^Ir but a;„ is close to 
Then, within a narrow band, |w„ — Qpt\ ~ Em, the Rabi 
oscillations develop an envelope whose shape is very sen- 
sitive to the detuning of the B^ driving frequency from 
Az- Experimentally, we have verified qualitatively and 
quantitatively these non-trivial modifications to Rabi nu- 
tation across a broad range of modulation-field frequency 
and amplitude with fairly simple NMR experiments on 
protons in water. 

When both and are smaller than rtj^, the Rabi 
oscillations are, in general, unaffected by the modulation, 
but even in this domain the oscillations acquire a fully 
developed envelope in the vicinities of certain modula- 
tion frequencies = 1)^/(2^ -I- 1). With increasing p 
the period of the envelope increases, while the domain, 
where it develops progressively shrinks, see Eq. 
( 45 ) . Physically, the limit p ^ 1 corresponds to a mul- 



tiphoton process where the role of photons is played by 
the "quanta of modulation" . This limit is easily captured 
within the adiabatic approximation^^. 

In certain semiconductor materials of interest, where 
magnetic resonance is detected electrically with the help 
of the pulsed techniqu^^H^ the effect of modulation can 
be even more peculiar. This is because the change of 
conductivity, measured in experiment, is determined by 
synchronized Rabi oscillations in pairs of spins, and re- 
flects the beating of these oscillations between the pair 
component^^. 

As a flnal remark, note that there is a conceptual sim- 
ilarity between the second-order Rabi oscillations, which 
develop at w„ = and Rabi-vibronic resonanc^^. In 
the latter case, the modulation of spacing between the 
Zeeman levels is accomplished via coupling of two-level 
system to a harmonic oscillator with frequency close to 



T he m ain result of the present paper is the mapping 
Eq. (17) which allows one to establish how a spin-i sys- 
tem under a weak resonant ac drive responds to a weak 
modulation of the longitudinal fleld with magnitude e„ 
and frequency w„ both smaller than the Zeeman split- 
ting Az- This response is strong and leads to a dramatic 
modification of the Rabi oscillations when £„, remain- 
ing much smaller than Az, exceeds the Rabi frequency, 
flu. Another instance when the response to modulation 
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